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Abstract 

Using Monte Carlo simulations, we study the properties of an elastic trian- 
gular lattice subject to a random background potential. As the cooling rate 
is reduced, we observe a rather sudden crossover between two different glass 
phases, one with exponential decay of translational correlations, the other 
with power-law decay. Contrary to predictions derived from continuum mod- 
els, no evidence of a crossover in the mean-square displacement, B(r), from 
quadratic growth at small r, to logarithmic growth at large r is found. 
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The structure of an elastic lattice in a random background has received much attention 
in recent years, due to its relevance to such systems as a Wigner crystal in a semiconductor 
with impurities ffl], a charge density wave in a weakly disordered medium 0, an atomic 
monolayer on an imperfect crystal surface ||, a magnetic bubble lattice in a ferromagnetic 
film with defects ||, and a vortex lattice in a disordered superconductor |J. The latter prob- 
lem has been investigated most intensively, both theoretically and experimentally, because 
translational correlations in the vortex lattice are responsible for such important properties 
of superconductors as resistivity and critical current. 

Most theoretical works on elastic lattices in a random background have concentrated on a 
continuum approach in which the deformation of the lattice is described by the displacement 
field u(r). In this model the background is described by a random potential V[r, u(r)] that 
satisfies V[r, u(r) + a»] = V[r, u(r)] which accounts for the periodicity of the lattice, (a* is 
a lattice vector.) The energy of this system is 

U = J d d r [dVuVu + V(t, u)] , (1) 

where the first term represents elasticity; aikim being the tensor of elastic moduli. When V = 
0, that is in the absence of pinning, the energy is minimized by u = const, which represents 
perfect translational order. In the presence of pinning the lattice develops deformations, 
u(r), and the question arises whether (and how rapidly) these deformations destroy long 
range translational correlations. A convenient measure of the disorder is 

B(r) =< [u(r) - u(0)] 2 > (2) 

If pinning is weak and only a small area of the lattice is of interest, the deformation u 
is small compared to the lattice spacing a. Then V(r, u) can be written as — f (r)-u and the 
problem reduces to the random force problem, for which the solution is known: a random 
force, no matter how weak, destroys the long range translational order in less than four 
dimensions. Simple statistical arguments yield ||[7j £>~(/ 2 /a 2 )r 4 ~ d , where f 2 denotes the 
variance of the random force. That is, Boer in three dimensions and Boer 2 in two dimensions. 
Equating B to a yields an estimate of the translational correlation length, £ oc f 2 l( d ~ 4 \ This 
estimate seems to be quite robust with respect to approximations and assumptions about 
the random potential. It has been pointed out, however, that the random field model cannot 
provide the correct rate at which translational correlations are destroyed for r > £, as it 
does not take into account the periodicity of the lattice. More sophisticated approaches 
based upon the Gaussian variational method || and the functional renormalization group 
H suggest that B oc (4 — d) In r for r > £ [|10| -[T2||. Another renormalization group approach 



|i~3|l suggests that the disordered lattice freezes below a certain temperature into a glass state 
with B oc In 2 r. The upshot of these analytical results is a much slower decay of translational 
correlations than predicted by the random field model 

The purpose of this paper is to compare the above predictions of the continuous model 
with Monte Carlo simulations of a discrete triangular lattice subject to a random background 
potential. As in the continuous approaches mentioned above, we limit our consideration to 
lattices free of dislocations. The results obtained under this assumption may be relevant 
to real systems since patterns of vortex lattices observed in decoration experiments show 
remarkably large areas free of dislocations. 
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We consider a two-dimensional triangular lattice of particles coupled by a harmonic, 
nearest-neighbor interaction, and subject to a static, random potential V(x, y). Let Xj = 
(xi, Hi) denote the position of particle i. Then the potential energy of the system is 

E = \ E(^-1) 2 + Enx,), , (3) 

<i,j> i 

where the first sum is over all nearest-neighbor pairs in the triangular lattice, and r^j = 
|x ; — Xj|. (All quantities are dimensionless in our formulation, with the basic length scale 
a = 1 set by the nearest-neighbor (NN) separation in the unstrained lattice and the basic 
energy scale e = 1/2 set by the energy of a NN pair with \r — 1| = 1, so that the spring 
constants have unit magnitude.) 

To avoid severe distortions of the lattice we impose a planarity constraint, which prevents 
any particle from escaping the "cage" defined by the current positions of its six nearest 
neighbors. (This is enforced by demanding that if x, y and z form a unit triangle in the 
unstrained lattice, then the angle between x-y and z-y not exceed ir.) We adopted this 
constraint in preference to equipping the particles with hard cores, since hard cores of a 
size sufficient to prevent this kind of distortion yield a rather strong anharmonicity. In the 
present model the anharmonicity is weak, with the deviation from equipartition amounting 
to < 5% at the temperatures of interest. 

The background potential V(x, y) is generated in two steps. We first generate R(i,j), 
an N x iV array (we used N = 101) of uncorrelated random numbers, uniformly distributed 
on [-1/2, 1/2]. We smooth this array by replacing each entry R(i,j) by the sum of the 
entry and its four nearest-neighbors, using periodic boundaries at the edges. The smoothing 
process is applied a total of three times, so that elements up to six units apart have a nonzero 
correlation. (Each entry of R is now a weighted sum of 25 random numbers, and so represents 
a good approximation to a Gaussian random variable.) To find V(x, y), we define a random, 
piece-wise constant (on a scale of ~ 10~ 3 ) map between points (x, y) and a set of four entries 
Ri, ...,i?4, where Rk = R(ik,jk)- The first entry is given by the modular expression: i\ = 
[7x (mod 101)], where the brackets denote the largest integer, and similarly for ji. Then %i = 
[fii (mod 100)], where / = 2049|i?i| is a random multiplier; j'2 is defined similarly, yielding 
i?2- We repeat the process, this time using / = 513|i?i + i^l, and then once more, with 
/ = 257|i?i + i?2 + -R3I - The resulting V(x,y) = A^ti R% is an (approximately) Gaussian 
random field with short-range correlations; C v (r) = (V(x)V(y))\ x - y \ =r /cr 2 drops from unity 
to about 0.3 for r ~ 10~ 3 , and then decays in roughly linear fashion, remaining essentially 
zero for r > 0.6. (We found that choosing prime factors in the multipliers hastened the decay 
of correlations.) The strength of the background potential is characterized by its standard 
deviation, a, which we control by varying the constant A. (a = 0.2 and 0.5 are used in 
this study. All simulations were performed on DEC Alpha workstations, and employed the 
random number generator supplied with the machine.) 

Thus the background potential differs from the kind typically employed in studies of 
off-lattice systems subject to quenched randomness which distributes a certain density of 
identical centers of force, with potential v(r), at random positions, q«, and sets V(x) = 
X)ii>(|x — Oil)- Our potential is a closer approximation to the Gaussian random field used 
in theoretical analyses. 

We simulated hexagonal-shaped lattices of M particles to a side (a total of 3M(M— 1) + 1 
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particles), with open boundaries. We report results for M = 60 and M = 120. We used 
open boundaries to eliminate global periodicity as a restraint on the growth of particle 
displacements. In each step of the simulation, a particle is selected at random and subjected 
to a trial displacement uniform on a square of side D = 0.5, symmetric about the origin. 
The move is accepted if the total change in energy AE < 0; if AE is positive the new 
position is accepted with probability e~ AE ^ T . Our time unit comprises one attempted move 
per particle. A preliminary study of the lattice without the random background potential 
(see also Ref. |fl4] ) revealed that the correlation function 

g G (r) = < e * G '(^)) , (4) 

(the thermal average is over all pairs with Xj — Xj = r in the unstrained lattice; G is a 
reciprocal lattice vector), shows a power-law decay, <?G( r ) ~ r ~ v with 77 proportional to 
temperature, T, as expected fllSf . 

Our primary interest is in the behavior of the mean-square displacement from equilib- 
rium, 

B(r) =< [x, - x, - if > , (5) 

where we average over all pairs with Xj — Xj - = r in the unstrained lattice, and the angle 
brackets denote an average over disorder. To avoid the effects of strong distortions that 
may appear at the boundary, we only average over particles at least M/2 sites distant from 
it. Configurations are generated by taking the perfect, unstrained lattice and permitting it 
to relax, in the presence of V(x), for = 10 4 — 10 5 timesteps at temperature 1. After 
this "annealing" phase the system is gradually cooled: the inverse temperature f3 = 1/T 
increases at a constant cooling rate T, until T = 0.01, at which time we compute B and 
other properties. (At this point the system is essentially at temperature zero, since the 
typical background energy a ^> T.) 

Depending on the cooling rate, we observe two qualitatively distinct kinds of B(r). In 
Fig. 1, for example, we show B(r) for M = 60, a = 0.2, = 10 4 , and a variety of different 
cooling rates. For relatively rapid cooling, (r > 4 x 10 -4 for the parameters of Fig. 1), B 
grows ~ lnr for small r before crossing over to some slower growth whose precise form is 
unclear. The correlation length £, defined via -B(£) = 1, decreases with T and in many cases 
exceeds the system size M. The logarithmic growth is reminiscent of a lattice in thermal 
equilibrium |T5|JT4| , suggesting that in this case the system has undergone a kind of glass 



transition, the final B{r) being a remnant of the thermal disorder when the lattice fell out 
of equilibrium with the heat bath. In this situation the background potential serves to 
"freeze" thermal fluctuations, even though the system does not have time to optimize its 
configuration with respect to V(r). Fig. 2 illustrates the similarity between B(r) in thermal 
equilibrium and in a rapidly cooled system. 

For lower cooling rates, the particles have the opportunity to explore more of the local 
potential energy landscape, and the lattice distortion, while more modest for small r, be- 
comes sizeable on large scales. In this regime the correlation length decreases as we reduce T; 
increased exploration results in greater distortion. An important open question concerns the 
nature of B(r) in the ground state: Does it grow ~ r 2 indefinitely, or cross over to a slower, 
perhaps logarithmic growth law? As we reduce the cooling rate, and hence probe nearer the 
ground state of the model, we find that B{r) maintains a faster than linear growth. Fig. 3 
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shows that in this regime, B(r) ~ r 2 to a good approximation. In no instance do we observe 
a crossover to logarithmic growth at large r, even though our data extend to r/£ ~ 14 for 
some systems. 

Fig. 4 shows the elastic and potential energies (u and v, respectively), per particle as 
a function of T, for the conditions of Fig. 1. While there are some fluctuations (due to 
computer-time restrictions we average over sets of only five independent disorder configu- 
rations in this series of studies), the energy decreases with cooling rate, as expected. We 
also plot B(r — 50) as a measure of the character of the final state (quenched thermal 
fluctuations versus optimized in the random background). It appears that the change in the 
behavior of B occurs rather suddenly, at T ~ 4 x 1CT 4 . 

In summary, we studied the mean-square displacement in an elastic lattice subject to 
a static random potential, and find no evidence of crossover to logarithmic growth in the 
well-relaxed regime. A crossover between two different glass phases, one with exponential 
decay of translational correlations, the other having power-law decay, is observed as the 
cooling rate is reduced. 

This work was supported by the Department of Energy under Grant No. DE-FG02- 
93ER45487. 
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Figure Captions 



FIG. I. B(r) versus r for M = 60, a = 0.2, = 10 4 , and various cooling rates V. 

FIG. 2. B(r) versus lnr for M = 60 and a = 0.2. Solid line: system in equilibrium at 
T = 0.3; broken line: T = 0.2; circles: system cooled to T = 0.01 at T = 1; squares: 
r = 0.5. 

FIG. 3. Mean-square displacement B(r) versus f = r/£. Dashed line: M = 120, a = 0.5, 
r = 5 x 10^ 4 ; solid line: M = 120, a = 0.2, T = 10" 3 ; dotted line: M = 60, a = 0.2, 
r = 10~ 4 . 

FIG. 4. Per-particle potential energy, v, elastic energy, u, and B(r = 50) versus V for the 
same parameters as in Fig. 1. 
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